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The nonlinear propagation of dust-acoustic (DA) waves in a magnetized dusty plasma with a pair 
of trapped ions is investigated. Starting from a set of hydrodynamic equations for massive dust fluids 
as well as kinetic Vlasov equations for ions, and applying the reductive perturbation technique, a 
Korteweg-de Vries (KdV)-like equation with a complex coefficient of nonlinearity is derived, which 
governs the evolution of small-amplitude DA waves in plasmas. The complex coefficient arises due 
to vortex-like distributions of both positive and negative ions. An analytical as well as numerical 
solution of the KdV equation are obtained and analyzed with the effects of external magnetic field, 
the dust pressure as well as different mass and temperatures of positive and negative ions. 


I. INTRODUCTION 

Recently, there has been a renewed interest in inves¬ 
tigating electrostatic disturbances in pair-plasmas and, 
in particular, plasmas with a pair of ions m- How¬ 
ever, nonthermal pair plasmas may frequently occur not 
only in semiconductors in the form of electron and ion 
holes, but also in many astrophysical environments, e.g., 
pulsars, magnetars, as well as in the early universe, ac¬ 
tive galactic nuclei and supernova remnants in the form 
of electrons and positrons [HHII]. On the other hand, 
a number of experiments have been conducted to cre¬ 
ate pair-ion plasmas using fullerene as ion source m- 
Furthermore, it has been observed that the dust parti¬ 
cles injected into a pair-ion plasma (e.g., /SFq plas¬ 
mas) can become positively charged when the number 
density of negative ions greatly exceeds that of electrons 
(tino ^ 500neo) uni HI- These pose some possibilities to 
investigate collective behaviors as well as the formation of 
nonlinear coherent structures in pair-ion plasmas under 
controlled conditions. The formation of phase space holes 
in pure pair-ion plasmas [S] as well as ion holes in dusty 
pair-ion plasmas |S] in the propagation of large amplitude 
electrostatic waves have been investigated in the recent 
past in which ions have been treated as trapped in self- 
created localized electrostatic potentials as prescribed by 
Schamel [15]. 

In this paper we present a theoretical study on the 
formation and the dynamics of small-amplitude solitary 
structures in a dusty plasma composed of charged dust 
particles and a pair of ions without electrons. In our the¬ 
oretical model the massive charged dusts are described 
by a set of fluid equations, while the dynamics of both 
positive and negative ions are governed by kinetic Vlasov 
equations. Using the reductive perturbation technique 
we show that the evolution of small-amplitude electro¬ 
static waves can be described by a Korteweg-de Vries 
(KdV)-like equation with a complex coefficient of the 


nonlinearity. A stationary as well as numerical solutions 
of the KdV equation are obtained and analyzed with the 
effects of external magnetic field, the dust pressure as 
well as different mass and temperatures of ions. 


II. BASIC EQUATIONS 

We consider the nonlinear propagation of dust-acoustic 
(DA) solitary waves in a magnetized dusty plasma which 
consists of positively or negatively charged mobile dusts 
and a pair of trapped ions with vortex-like distributions. 
The dust particles are assumed to have equal mass and 
constant charge. The collisions of all particles are also 
neglected compared to the dust plasma period. Fur¬ 
thermore, in dusty pair-ion plasmas the ratio of electric 
charge to mass of dust particles remains much smaller 
than those of positive and negative ions. We also assume 
that the size of the dust grains is small compared to the 
average interparticle distance. The static magnetic field 
is considered along the z-axis, i.e., B = Bqz. While the 
dynamics of massive charged dusts in the propagation of 
DA waves (vtd ^ Vp ^ Vp^n, where Vtj (= ^ksTj/rrij) is 
the thermal velocity of j-th species particles and Vp is the 
phase velocity of the wave) is described by a set of fluid 
equations ([^ and ([^, the dynamics of singly charged 
positive and negative ions are described by the Vlasov 
equations 

^-f V • (ridVd) = 0, (1) 


^ + (vd • V) Vd = ^ (E + Vd X Bo) 
at nid 


VP 
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The system of equations is then closed by the Poisson 
equation 
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In Eqs. 0-0, Qj, rij, Vj, fj and rrij respectively, de¬ 
note the charge, number density (with its equilibrium 
value rijo), velocity, velocity distribution function, and 
mass of j-species particles. Also, qd = azde with a = ± 
denoting for positively/negatively charged dusts and Zd 
the charge state. Also, E = — V(/) is the electric field with 
(j) denoting the electrostatic potential and P is the dust 
thermal pressure given by the adiabatic equation of state 
P/Po = {nd/udo)^■ Here, 7 = 5/3 is the adiabatic index 
for three-dimensional configuration and Pq = ridoksTd is 
the equilibrium dust pressure with kg denoting the Boltz¬ 
mann constant and Tj the thermodynamic temperature 
of j-species particles. Furthermore, the ion densities are 
given by 



In what follows, we recast Eqs. 0-0 in terms 
of dimensionless variables. To this end the physical 
quantities are normalized according to rij —> 

(v, Vd) (v, 'Vd)lcd, 4‘ ecji/ksTp with e denot¬ 
ing the elementary charge, fj —)■ fjVtp/rijo, where Cd = 
yjZdkBTp/md = oJpd^D is the dust-acoustic speed with 
ujpd = yj-iTTUdoz^e^/rud and \d = y/ksTp/inudoZde'^ de¬ 
noting, respectively, the dust plasma frequency and the 
plasma Debye length. The space and time variables are 
normalized by Ad and respectively. Thus, from Eqs. 
0-0 we have following set of normalized equations. 

^ -b V • (udVd) = 0, ( 6 ) 

+ aV(l) = auJc-Vd X z - ^adu'^^^^Vud, (7) 


are constant of motion of the Vlasov Eq. ([^, are cho¬ 
sen m for the excitation of localized solitary waves so 
that (i) they are continuous, and both the free particle 
distributions are Maxwellian distribution where ^ 0 

at 1^1 —>■ ±00 and trapped particles are absent, (ii) both 
trapped particle distributions are Maxwellian (with also 
negative temperatures). Thus, fj (for free and trapped 
particles) are (with a suitable choice of the normalization 
constants) [TKHTI^ for positive ions 


fpfiv) 


1 
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and for negative ions 
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exp 
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\v\ < \/2ip/r 

(15) 

where m (= nin/mp > 1 ) is the mass ratio, cr (= Tp/Tn > 
1 ) is the temperature ratio and aj, for j = p, n, measure 
the inverse of the trapped positive and negative ion tem¬ 
peratures which may be negative {aj < 0 ) correspond¬ 
ing to a depression in the trapped particle distribution. 
The case of aj —>■ 0 represents the plateau (constant or 
flat-topped) and aj —)■ 1 corresponds to the Boltzmann 
distribution of ions. Next, integrating the particle dis¬ 


tribution functions (12)-(15) over the velocity space, i.e., 
using Eq. (10) we obtain the number densities rij for 


positive and negative ions as 


(5 


dfj 

dt 
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( 10 ) 


np{(j)) = /(-(/)-b 



e erf {yj-aptf) , ap > 0 

ap(j)^ , CTp < 0, 

(16) 


> 0 

< 0 , 

(17) 

where I{x) = exp(a;) [1 — erf(y^)]. The error and Daw¬ 
son functions erf(a;) and W(x) are, respectively, given by 


nn{4>) = H^4>) 



erf (v^ 


</) , 


where d/dt = d/dt -b • V, a = ±1 for posi¬ 
tively/negatively charged dusts, Wc = IqdlBo/mdtOpd is 
the dust-cyclotron frequency normalized by the dust 
plasma frequency, ad = Td/TpZd, 5 = y/z^m^Jmd, 
Q = ±1 for positive/negative ions and = rijo/Zdno 
are the density ratios (j = p, n) which satisfy the follow¬ 
ing charge neutrality condition at equilibrium: 

Pp -bo P'n- (11) 


2^2 2^2 
erf(a;) = —j= / e * dt, W{x) = e ^ / e* dt. (18) 
Jo Jo 

In the small amplitude limit 0 <C 1 , so that acj) 1, we 
obtain from Eqs. ([T^ and jEL the following expressions 
for the number densities nsHTi] 


..^1 J. 4(1 - ap) -iji 

rip « 1 - 0 - „ ^ (- 0 ) 

oy TT 


(19) 


We neglect the ion inertial effects compared to the 
charged dusts, i.e., <5 —)■ 0 in Eq. The distribu¬ 

tion functions fj for positive and negative ions, which 


« 1 -b (cr0) - (o-0)^^^ + ^ (cr0)^ . (20) 
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FIG. 1. Profiles of |$| given by Eq. ( |33[ ) are shown with respect to ^ (and with cr ~ 1) for different values of the plasma 
parameters as in the figure. The Hxed parameter values for the subplots (a) to (d), respectively, are (a) a = 1, jj,p = 0.2, 1^ = 
0.1, uic = 0.5, ad = 0.06 and uo = 0.1, (b) a = 1, ^lp — 0.2, G = 0.1, cTp = cr„ = 0.3, ad = 0.06 and uq = 0.1, (c) q = 1, Up = 
0.2, Iz = 0.1, ap = an = 0.3, uJc = 0.5, ad = 0.06 and uo = 0.1, and (d) a = 1, ^.p = 0.2, ap = an = 0.3, uJc = 0.5, ad = 0.06 
and Mo = 0.1. 



FIG. 2. The space-time evolution of the soliton profile |‘I>| [numerical solution of Eq. @] is shown at r = 100 [Subplots (a) 
and (b)] and r = 160 [Subplots (c) and (d)|. While the plots (a) and (c) show the solton profiles with space, plots (b) and are 
the corresponding contour plots. The parameter values are the same as for the dashed line in Eig. m and (7 ~ 1. 


III. EVOLUTION EQUATION where e is a small parameter measuring the strength of 

nonlinearity. The dependent variables are expanded as 

In order to derive the evolution equation for the DA 
waves, we transform the space and time variables accord¬ 
ing to [T^ 


^ {IxX + lyU + IzZ - Mt ), T = (21) 
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n = 1 + H— ■ 

^ = e^d) + e3/2^(2) ^ .., 


■u^ = 1 + + e' 


‘^a:.47 — ^ 




,3/2^(2) 


( 22 ) 


The anisotropy in Eq. (22) for the transverse velocity 


components of dust fluids is introduced on the assump¬ 
tion that the dust gyromotion is a higher-order effect 
than the motion parallel to the magnetic field. Next, 
we substitute Eqs. (21) and (22) into Eqs. ([^-([^ and 
equate different powers of e successively. In the lowest 
order we obtain the following first-order quanti¬ 

ties 


(23) 


,( 1 ) = 


di) = 

x,y 


CUc 




(24) 


and the dispersion relation for the nonlinear wave speed 
given by 


M = l, 


5 1 

H-;- 

3 iJ.p + ann 


1/2 


(25) 


Replacing M by w/fc, one can obtain the same dispersion 
relation after Eourier analyzing the linearized basic equa¬ 
tions (§-(§, i.e., assuming the perturbations as oscilla¬ 
tions with the wave frequency w and the wave number 
k. We find that the phase speed M (normalized by the 
DIA speed Cd) can be larger or smaller than the unity de¬ 
pending on the choice of the parameter values. The value 
of M increases with increasing values of both and (Td- 
However, its values can slowly decrease with increasing 
values of the density ratios fij as well as the temperature 
ratio a. From Eq. ([^, collecting the coefficients of 
we obtain 


27| a nd (28), and the coefficient of vanishes 
by Eq! (|25[). Thus, arranging the terms and using Eq. 


(26), 


(23) one obtains the following KdV-like equation 


54> 


+ ^Ap\/—4* + + B — 0, (29) 

where d* = It follows that Eq. (29) has a complex 


^3$ 


solution for $. Typically, for <I> ~ rexp(z0), where r is 
a real function of ^ and r, and 0 is a constant, one can 
have = i'/^. Thus, Eq. (29) can be written as 


^ + A^- + B^ = 0 . 


(30) 


where the coefficients of nonlinearity {A = An + iAp) and 
dispersion (i?) are given by 


A, = 


(l-crj)/ij /Tj 


M{pp + a^n) \T, 


3/2 


(31) 


B = 


2M 


11 

a,?rt 


CTfin 


(32) 


The nonlinear coefficient A becomes complex due to 
vortex-like distributions of two oppositely charged par¬ 
ticles. In absence of one of them, A becomes real, and 
one can then obtain solitary waves with positive or neg¬ 
ative potential. A stationary soliton solution of Eq. (30) 
can easily be obtained with its absolute value as (Fbr 
details see Appendix A) 


|$| = $0 sech'‘ [(^ - uqt) jW ], 


(33) 


where uq is a constant, and $o = (15mo/8|A|)^ and W = 
a/16H/mo are the amplitude and width of the soliton 
respectively. 


IV. RESULTS AND DISCUSSION 


M 






dr 


Y. h 




(26) 


Similarly, equating the coefficients of from the x- 
and y-components of Eq. ([^, and the coefficients of 
from the z-component of Eq. 0 we successively obtain 


av 


( 2 ) 

x,y 


= ± 


AI dv'^yl 

oJc 


Ml 




UJ 
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+ o^d (Mp + cr^n) 


de 

(27) 


M 


a? 




(28) 


From the coefficients of of Eq. (|^, we obtain an 
equation in which is eliminated by the use of Eqs. 


We numerically analyze the solution ( |33[ ) with dif¬ 
ferent plasma parameters as shown in Fig. Since 
(j = P: ^) represents the reciprocal temperature of the 
trapped positive and negative ions, and can be allowed 
from their negative to positive values corresponding to 
different trapped particle distributions, we consider neg¬ 
ative, zero as well as positive values of aj. 

From Fig. [^a), it is seen that as cxj increases from 
aj = 0 (corresponding to a constant or flat-topped dis¬ 
tribution of ions) to aj ~ I (corresponding to the Boltz¬ 
mann distributions of ions), both the amplitude and 
width of the soliton increase (See the solid and dashed 
lines). Note here that the values of aj > I, for which the 
influence of the trapped ions are inverted, may be physi¬ 
cally unrealistic as those correspond to a more steepened 
wave which can become unstable due to more peaked 
bump of the ion distributions. However, as the absolute 
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value of A starts increasing for aj < 0, which corresponds 
to a depression in the trapped particle distribution, both 
the amplitude and width of the soliton are reduced (See 
the dotted line). The same can further be enhanced for 
values of aj satisfying crpcr„ < 0 (See the dash-dotted 
line). 

Figure [^b) shows the soliton profile with the influ¬ 
ence of the external magnetic field. Since oJc contributes 
only to the dispersive coefhcient B of Eq. p0| ), the ef¬ 
fect of the magnetic field with increasing its intensity is 
to reduce the width (without changing the amplitude) 
of the soliton. Thus, the external magnetic field makes 
the solitary structure more spiky. However, for stronger 
magnetic fields with Wc ^ 1, the width remains almost 
unaltered as in this case B ~ l‘l/2M. 

The thermal effects of charged dusts are shown in Fig. 
[^c) . It is found that the effect of the dust thermal pres¬ 
sure (Td is to enhance both the amplitude and width of 
the soliton. The enhancement is due to the fact that as 
(Jd increases, the values of |H| {B) decrease (increase), 
and hence the increase in both the amplitude and width. 
However, an opposite trend occurs by the effects of the 
positive to negative ion temperature ratio a (not shown 
in the figure). Typically, it reduces both the soliton am¬ 
plitude and width significantly with a small increment of 
its value. 

Figure [^d) exhibits the effects of the obliqueness of 
propagation 4 and the relative (to dusts) concentration 
of positive ions /ip. We find that both the amplitude 
and width of the soliton are greatly enhanced by a small 
increment of Iz [Since Aj (B) is inversely (directly) pro¬ 
portional to Iz]. However, as the positive ion concentra¬ 
tion increases, the amplitude gets reduced but the width 
increases. 

Next, we numerically solve Eq. (30) by the Runge- 
Kutta scheme with an initial condition of the form d)(^) = 
0.05 sech‘^(f/10) exp(—j^/15) and time step dr — 0.001. 
The development of the wave form |$| after a finite in¬ 
terval of time is shown in Fig. The parameter values 
are considered as the same as for the dashed line in Fig. 
m- It is seen that the leading part of the initial wave 
steepens due to positive nonlinearity. As the time goes 
on the pulse separates into solitons and a residue due to 
the wave dispersion [See the subplots (a) and (b)]. It is 
found that once the solitons are formed and separated, 
they propagate in the forward direction without chang¬ 
ing their shape due to the nice balance of the nonlinearity 
and dispersion [See the subplots (c) and (d)]. 


V. CONCLUSION 

We have investigated the nonlinear propagation of 
dust-acoustic waves in a magnetized plasma which con¬ 
sists of warm positively charged dusts and a pair of free, 
as well as, trapped ions. We have shown that the evo¬ 
lution of small-amplitude DA waves can be described by 
a KdV-type equation with a complex coefficient of the 


nonlinearity. Such complex coefhcient appears due to 
vortex-like distributions of both the ion species. The 
KdV equation is solved both analytically and numeri¬ 
cally. The properties of the absolute value of $ are only 
exhibited graphically. It is shown that while the external 
magnetic held only inhuences the width of the soliton, 
the trapped ion temperatures, the thermal pressures of 
ions and dusts, the relative concentration of positive ions 
as well as the obliqueness of propagation have signihcant 
effects on both the amplitude and width of the solitons. 
We stress that other solutions [TM2^ than those pre¬ 
sented here of the complex KdV equation could of interest 
but beyond the scope of the present work. To conclude, 
the present results should be useful in understanding the 
nonlinear features of electrostatic localized disturbances 
in laboratory and space plasmas. 


ACKNOWLEDGMENTS 

This work was partially supported by the SAP-DRS 
(Phase-H), UGC, New Delhi, through sanction letter No. 
F.510/4/DRS/2009 (SAP-I) dated 13 Oct., 2009, and 
by the Visva-Bharati University, Santiniketan-731 235, 
through Memo No. REG/Notice/156 dated 07 January, 
2014. APM thanks Dr. M. M. Panja of Department 
of Mathematics, Visva-Bharati, Santiniketan, India for 
some useful discussions. 


Appendix A: Stationary solution of the KdV-like 
equation 


Equation (301 is recast as 




(Al) 


Next, we apply the transformation rj = ^ — uqt to obtain 
from Eq. (Al) 


drj 


— ( 5$ - uo$ -I- - 1=0 




(A2) 


where the dot denotes differentiation with respect to rj. 
Integrating Eq. (A2) with respect to rj and using the 
boundary conditions 4>, $ —>• 0 as ^ ±00 we get 


- Mo® -I- = 0. 

3 


Multiplying Eq. 
respect to r/, we obtain 


(A3) 

by 2 <i> and integrating once with 


= 0 , 

15 


(A4) 


where we have used the boundary conditions $, d* —>■ 0 . 
From Eq. (A4) we have 




(A5) 
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or, 


’ / /- 

•’ (8A/15B)y$ 

which gives (a = mo/-S and b = %A/lbB) 


_ y (A6) Thus, we obtain a soliton solution of Eq. (30) as 


l_tanh2( 


(A9) 






(A7) 


or, 



=p tanh 



(A8) 


or, $ = 



2 

sech^ 



(AlO) 
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